function [ y,dy ] = sigmoid( x )
  y = sigmf(x,[1,0]);
  dy = y.*(1-y);
end  % sigmoid
